Charged-particle decay and suppression of small-scale power 
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We study the suppression of the small-scale power spectrum due to the decay of charged matter 
to dark matter prior to recombination. Prior to decay, the charged particles couple to the photon- 
baryon fluid and participate in its acoustic oscillations. However, after decaying to neutral dark 
matter the photon-baryon fluid is coupled only gravitationally to the newly-created dark matter. 
This generically leads to suppression of power on length scales that enter the horizon prior to 
decay. For decay times of ~3.5 years this leads to suppression of power on subgalactic scales, 
bringing the observed number of Galactic substructures in line with observation. Decay times 
of a few years are possible if the dark matter is purely gravitationally interacting, such as the 
gravitino in supersymmetric models or a massive Kaluza-Klein graviton in models with universal 
extra dimensions. 
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The standard inflation-inspired cosmological model, 
with its nearly scale-invariant power spectrum of primor- 
dial perturbations, is in remarkable agreement with ob- 
servation. It predicts correctly the detailed pattern of 
temperature anisotropics in the cosmic microwave back- 
ground (CMB) 0, and accurately describes the large 
scale clustering of matter in the Universe 0. However, 
on subgalactic scales there are possible problems with the 
standard cosmology that warrant further investigation. 
Namely, the model overpredicts the number of subgalac- 
tic halos by an order of magnitude compared to the 11 
observed dwarf satellite galaxies of the Milky Way 0. 
Several possible resolutions have been proposed to this 
apparent discrepancy, ranging from astrophysical mech- 
anisms that suppress dwarf-galaxy formation in subgalac- 
tic halos (see, for example, Ref. jj]) to features in the in- 
flaton potential that suppress small-scale power and thus 
reduce the predicted number of subgalactic halos 

In this Letter, we show that if dark matter is produced 
by the out-of-equilibrium decay of a long-lived charged 
particle, then power will be suppressed on scales smaller 
than the horizon at the decay epoch. Unlike some other 
recent proposals, which suppress small-scale power by 
modifying the dark-matter particle properties 6], ours 
modifies the dark-matter production mechanism. In the 
model we discuss here, prior to decay, the charged parti- 
cles couple electromagnetically to the primordial plasma 
and participate in its acoustic oscillations. After decay, 
the photon-baryon fluid is coupled only gravitationally to 
the neutral dark matter. This generically leads to sup- 
pression of power for scales that enter the horizon prior 
to decay. This suppression, reduces the amount of halo 
substructure on galactic scales while preserving the suc- 
cesses of the standard hierarchical-clustering paradigm 
on larger scales. Apart from the changes to the model 
due to the decay process, we adopt the standard flat- 
geometry ACDM cosmological model with present-day 
dark-matter density (in units of the critical density) = 
0.25, baryon density VL\, = 0.05, cosmological constant 



tt A = 0.70, Hubble parameter H = 72 kms^Mpc^ 1 , 
and spectral index n = 1. 

In the standard ACDM model the initial curvature per- 
turbations of the Universe, presumably produced by in- 
flation or some inflation-like mechanism, are adiabatic 
(perturbations in the total density but not the relative 
density between species) and Gaussian with a nearly 
scale- invariant spectrum of amplitudes. These initial per- 
turbations grow and react under the influence of grav- 
ity and other forces, with the exact nature of their be- 
havior dependent upon the species in question. Because 
dark-matter particles are, by assumption, cold and colli- 
sionless the fractional dark-matter-density perturbation 
5d = Spd/pd can only grow under the influence of gravity. 
The baryonic species, being charged, are tightly coupled 
by Coulomb scattering to the electrons, which are them- 
selves tightly coupled to the photons via Thomson scat- 
tering. The baryons and photons can thus be described 
at early times as a single baryon-photon fluid, with the 
photons providing most of the pressure and inertia and 
the baryons providing only inertia. Gravity will tend to 
compress this baryon-photon fluid, while the radiation 
pressure will support it against this compression. The 
result is acoustic oscillations, and the baryon density per- 
turbation 8b = Spb/pt and photon density perturbation 
<5 7 = Spy/py will oscillate in time for length scales inside 
the horizon (on length scales larger than the horizon the 
pressure can have no effect). At early times these per- 
turbations are very small and linear perturbation theory 
can be applied. This allows an arbitrary density field 
to be decomposed into a set of independently evolving 
Fourier modes, labeled by a wavenumber k. Fig.^shows 
the growth of dark-matter perturbations under the influ- 
ence of gravity, and the oscillatory behavior of the baryon 
perturbation for the same wavenumber. 

We choose to work in the synchronous gauge where the 
time slicing is fixed to surfaces of constant proper time 
so that particle decays proceed everywhere at the same 
rate. In the synchronous gauge the standard linearized 
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FIG. 1: The evolution of the comoving wavenumber k = 
3.0 Mpc -1 density perturbations in the early Universe for 
dark matter (dashed line) and baryons (dotted line). The 
dark-matter perturbation always grows under the influence 
of gravity while the baryonic perturbation oscillates due to a 
competition between gravity and the photon pressure. 



linearized equations. The last terms on the right-hand- 
sides of Eqs. and (J5J account for Thomson scattering 
between baryons and photons, and are responsible for 
keeping them tightly coupled in the early Universe. In 
these equations ot is the Thomson cross section, n e is 
the electron number density, and c s is the intrinsic sound 
speed of the baryons. During tight coupling the second 
moment 7 of the photon distribution and other higher 
moments can be neglected, and the radiation can reliably 
be given the fluid description described above. 

We now show how Eqs. are modified by the 

decay of a long-lived metastable charged particle to dark 
matter in the early Universe. We assume that the decay 
is of the form g ± — ► i^d, so the decay of each charged par- 
ticle q ± produces a dark-matter particle d and a charged 
lcpton ^ . Denoting the decaying charged component 
by the subscript ( q\ the background density p q evolves 
according to the equation, 



Pq 



a a 

-3-pq Pq , 

a t 



(6) 



where r is the lifetime of q . The first term just accounts 
for the normal a~ 3 scaling of non-relativistic matter in 
an expanding universe, while the second leads to the ex- 
pected exponential decay of the comoving density. For 
the dark matter we have 



a a 
Pd = -3-Pd + \-p q , 



(7) 



evolution equations for perturbations in Fourier space are 
(e.g., 0) 



where A = m c i/m q is the ratio of the masses of dark 
matter particle the charged particle. The energy density 
in photons evolves according to 



Sd — — 9d — r ^ , 9d — ——9d j 
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and 



9 7 = k ( —S-y — 7 ) + an e <TT(6b — Q-y) 
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Pi = -4-p 7 + (1 - A)~Pq ■ 
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where Ob, Od, and # 7 are the divergence of the baryon, 
dark-matter, and photon fluid velocities respectively and 
an overdot represents a derivative with respect to the 
conformal time r\. Here h is the trace of the spatial met- 
ric perturbations hij. Its evolution is described by the 
linearized Einstein equations, which close this system of 



This last equation follows from the assumption that 
the produced lepton initiates an electromagnetic cascade 
which rapidly (compared to the expansion timescale) 
thermalizes with the photon distribution. In practice the 
last term on the right-hand-side of Eq. JSJ is negligibly 
small because the decay takes place during the radiation 
dominated era when p 1 ^> p q . Furthermore, limits on the 
magnitude of ^-distortions to the blackbody spectrum of 
the CMB constrain |1 — A| to be a small number, as we 
discuss below. 

Using covariant generalizations of Eqs. we can 

derive how Eqs. are modified by the transfer 

of energy and momentum from the component to 
the dark matter during the decay process. Since the 
charged '<?' component and the baryons are tightly cou- 
pled via Coulomb scattering they share a common veloc- 
ity Op = Ob = q . This makes it useful to describe them in 
terms of a total charged-species component with energy 
density pp = p b + p q , which we denote here by the sub- 
script '/?'. Because in the synchronous gauge the decay 
proceeds everywhere at the same rate this description is 
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FIG. 2: The evolution of the comoving wavenumber k — 30.0 Mpc -1 (left panel), k — 3.0 Mpc -1 (center panel), and 
k — 0.3 Mpc -1 density perturbations in the early Universe for dark matter in the ACDM model (dashed line) and in the model 
with r = 3.5 yr (solid line). The '/?' perturbation is represented by the dotted line. Due to being sourced by the low amplitude 
'/?' perturbations at early times the dark matter perturbation in the model with a decaying charged component is suppressed 
relative to the standard ACDM case for k = 3.0 Mpc" 1 . For k ^> 3.0 Mpc -1 (very small scales) 8d tracks the oscillations in Sp 
before decay, while for k <C 3.0 Mpc -1 (large scales) 5d follows the standard growing evolution. 



even more useful as 5p = 5j, = 5 q is maintained at all 
times for adiabatic initial conditions. In terms of these 
'/?' variables, then, we have 

Sd = -0d - ~h + \—-{5p - 5 d ) , (9) 

2 PdT 

0d = --Od + X^-(9 fJ -O d ), (10) 
a PdT 

h = -h-\h, (11) 

Bp = --0/3 + c 2 k 2 5p + ^<m e a T (0 7 - 9p) , (12) 
a 6pp 

5, = -^8 1 -h+(l-\)?l-(5 -6 1 ), (13) 

O O P-y T 

and 

9*y = k 2 ( — <5 7 — 7 ^ + an e GT{0p — # 7 ) 

+ (1 _^(^_ fl7 ). (14) 

We now describe how small-scale modes that enter 
the horizon prior to decay are suppressed relative to 
those modes that enter the horizon after decay. Due 
to the Thomson collision terms the '/?' component and 



the photons will be tightly coupled as a '/^'-photon fluid 
at early times and this fluid will support acoustic oscil- 
lations. Furthermore, Eqs. © and QlUf) show that the 
dark-matter perturbations are strongly sourced by the 
perturbations of the '/3' component prior to decay, when 
the ratio p q /pd is large. Dark-matter modes that enter 
the horizon prior to decay will thus track the oscillations 
of the '/3'-photon fluid rather than simply growing un- 
der the influence of gravity. After decay, when the ratio 
Pq/Pd is small, the source term shuts off and dark-matter 
modes that enter the horizon undergo the standard grow- 
ing evolution. In Fig. [21 we follow the evolution of the 
dark-matter perturbations through the epoch of decay. 
We modified CMBFAST Q to carry out these calculations. 

In order to suppress power on subgalactic scales the 
decay lifetime must be roughly the age of the Universe 
when the mass enclosed in the Hubble volume is equal to 
a galaxy mass; this occurs when r ~ years. In Fig. [3] we 
plot the linear power spectrum of matter density fluctu- 
ations at the present day for a charged-particle lifetime 
r = 3.5 yr assuming a scale-invariant primordial power 
spectrum. We see that power is suppressed on scales 
smaller than ~ 0.3 Mpc relative to the standard 
ACDM power spectrum. Suppression of power on these 
length scales reduces the expected number of subgalactic 
halos, bringing the predictions in line with observation 
5] without violating constraints from the Lyman-alpha 
forest 0. Of course, the model reproduces the successes 
of the standard ACDM model on larger scales and in the 
CMB. 

The requirements of the charged-particle species are 
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FIG. 3: The linear-theory power spectrum of matter density 
fluctuations in the standard ACDM model (dashed line), and 
in the charged decay to dark matter model (solid line) with 
r = 3.5 yr. The charged decay model matches the standard 
ACDM model on length scales larger than 0.3 Mpc, but power 
drops of sharply below 0.3 Mpc. 

that it have a comoving mass density equal to the dark- 
matter density today and have a lifetime of r ~ 3.5 yr. 
In order to satisfy the constraint to the CMB chemi- 
cal potential [10(, the fractional mass difference between 
the charged and neutral particles must be Am/m < 
3.6 x 1CV 3 , and in order for the decay to be allowed 
kinematically the mass difference must be greater than 
the electron mass. One possibility is the SuperWIMP 
scenario of Ref. [TH ] in which a charged particle may 
decay to an exclusively gravitationally interacting par- 
ticle. For example, in supersymmetric models, thejdecay 
of a selectron to an electron and gravitino e — > e G with 
nie ~ nig > 122 TeV would satisfy these constraints, as 
would the decay of a KK-electron to an electron and KK- 
graviton e 1 — > e G 1 with m e i rs m G i > 72 TeV in the 
case of the single universal extra dimension Kaluza-Klein 
(KK) model discussed in Refs. [H03- 

Such masses are 

larg er than the unitary bound for thermal production 
|l3| . but might be accommodated through nonthermal 
mechanisms or if the next-to-lightest partner is a squark 
which might then interact more strongly and thus evade 
this bound. There may also be viable scenarios involving 
nearly-degenerate charged and neutral higgsinos. 

It should be noted that the recent WMAP evidence 
for early star formation 

El 

argues against the suppres- 
sion of small-scale power, but these results are not yet 



conclusive. If it does turn out that traditional astrophys- 
ical mechanisms can explain the dearth of dwarf galaxies, 
then our arguments can be turned around to provide con- 
straints to an otherwise inaccessible region of the param- 
eter space for decaying dark matter 14] . Finally, if the 
mechanism we propose here is realized in nature, then 
the dearth of small-scale power, along with the detection 
of a non-zero CMB chemical potential, would be a pow- 
erful probe of the particle spectrum of the new physics 
responsible for dark matter. 
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